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ABSTRACT 

Recent laboratory experiments have demonstrated that, even though contribution 
from other reaction channels cannot be neglected, unsaturated hydrocarbons easily 
break their multiple C-C bonds to form CO after their interactions with atomic oxy- 
gen. Here we present an upgraded chemical modelling including a revision of the 
reactions between oxygen atoms and small unsaturated hydrocarbons for different as- 
trochemical environments. A first conclusion is that towards hot cores/corinos atomic 
oxygen easily degrades unsaturated hydrocarbons directly to CO or to its precursor 
species (such as HCCO or HCO) and destroys the double or triple bond of alkenes 
and alkynes. Therefore, environments rich in atomic oxygen at a relatively high tem- 
perature are not expected to be rich in large unsaturated hydrocarbons or polycyclic 
aromatic hydrocarbons. On the contrary, in O-poor and C-rich objects, hydrocarbon 
growth can occur to a large extent. In addition, new radical species, namely ketyl and 
vinoxy radicals, generated from other reaction channels can influence the abundances 
of hydrocarbons towards hot cores. We, therefore, suggest they should be included in 
the available databases. Hydrocarbon column densities are calculated in the 10 13 -10 15 
cm~ 2 range, in good agreement with their observed values, despite the small number 
of data currently published in the literature. 
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1 INTRODUCTION 

The interstellar medium (ISM) is extremely heterogeneous 
with a broad range of temperatures, densities and extinc- 
tions (Snow & Bierbaum 2008). As a consequence a wide 
variety of chemical processes can occur. Oxygen is an impor- 
tant player in the ISM chemistry because it is the most abun- 
dant element after hydrogen and helium. Interstellar oxygen 
can be found in different forms and compounds, such as neu- 
tral or ionised gas, in the water of icy mantles or even in 
refractory materials. Its gas-phase abundances in the diffuse 
medium are well known (Meyer et al. 1997), while the oxygen 
chemistry in denser environments seems to be poorly under- 
stood due to observational constraints (Jensen et al. 2005). 
On the other hand oxygen chemistry has been widely inves- 
tigated towards extragalactic and galactic environments: in 
line with the fact that the C/O ratio determines the kind 
of reactions taking place in a specific region of ISM (Tsuji 
1973), Agundez & Cernicharo (2006), looking at the circum- 
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stellar shells of IRC +10216, observed that carbon is locked 
either as CO or other O-bearing species when C/O < 1, 
while a higher concentration of reduced C-bearing species is 
present when C/O > 1. The latter condition leads to an in- 
crease in the growth of unsaturated hydrocarbons. This has 
been extensively confirmed by modelling and observations 
and in particular Woods et al. (2012) found a large presence 
of unsaturated hydrocarbons, such as acetylene, by comput- 
ing the physical conditions around carbon-rich AGB stars 
at sub-solar metallicities. 

Unsaturated hydrocarbons are organic molecules char- 
acterised by double or triple covalent bonds between ad- 
jacent carbon atoms; presently, there is no generally ac- 
cepted explanation of their formation mechanisms but, as 
mentioned above, their abundances seem to be inversely cor- 
related to the availability of atomic oxygen (Fortney 2012). 
Their most likely reactions start with the addition of atomic 
and radical species to their multiple bonds. In most cases, 
the elimination of a hydrogen atom or of a group follows, 
leading to more complex compounds (see, for instance, the 
cases of the reactions C + C2H2, C2 + C2H2, CN + C2H2, 
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Kaiser 2000, Leonori et al. 2008a, Kaiser et al. 2003, Leonori 
et al. 2008b, Huang et al. 2000, Leonori et al. 2010). Recent 
experiments in the gas-phase (Capozza et al. 2004; Fu et al. 
2012a; Fu et al. 2012b; Leonori et al. 2012) suggest that the 
C-C bond breaking channel with the formation of CO is an 
important reaction channel for the reaction O + C2H2 and 
O + C2H4, while it is by far the dominant loss path for the 
reaction O + CH2CCH2 (allene). The fact that the reac- 
tions of atomic oxygen with unsaturated hydrocarbons lead 
in one step to CO can have a pivotal role in estimating the 
CO abundances both in the gas-phase and on grain surfaces. 
Carbon monoxide is in turn a key species for the synthesis 
of organic complex molecules on grain surfaces since it can 
easily hydrogenate to methanol as experimentally demon- 
strated by Modica & Palumbo (2010). Another important 
aspect regarding these small hydrocarbons is their loss paths 
to form CH2. The latter molecular radical is indeed known to 
be a potential precursor of polycyclic aromatic hydrocarbons 
(PAHs) and many studies have been performed to stress 
several implications due to the presence of these aromatic 
hydrocarbons on the dust (Verstraete 2011; Salama 2008); 
on the other hand small unsaturated hydrocarbons are the 
building blocks of polyynes (organic compounds with alter- 
nating C-C single and triple bonds) and cumulenes (hydro- 
carbons with three or more consecutive C-C double bonds) 
that are also known as starting points for the synthesis of 
PAHs. 

In the present study we focus on the reactions between 
O atoms and small unsaturated hydrocarbons, namely, 
acetylene, ethylene, methylacetylene and allene, under in- 
terstellar medium conditions. In particular, we propose a 
revised chemical scheme (compared to those included in the 
available databases) where our reaction rates are based on 
new experimental branching ratios (BRs) for the different re- 
action channels (Capozza et al. 2004; Fu et al. 2012a; Fu et 
al. 2012b; Leonori et al. 2012). In order to assess the effects of 
the recent laboratory data on the hydrocarbon abundances 
across different astrophysical environments, we model four 
typical astrochemical regions and we compare our theoreti- 
cal results with the observations. The paper is organised as 
follows: Section 2 describes the laboratory technique used 
and the revised chemistry we inserted in our model; Sec- 
tion 3 focuses on the modelling of four different astrochemi- 
cal environments; results are discussed in Section 4 and our 
conclusions are summarised in Section 5. 



2 THE REACTIONS OF ATOMIC OXYGEN 
WITH SMALL UNSATURATED 
HYDROCARBONS 

The following gas-phase neutral-neutral reactions have been 
experimentally investigated at the University of Perugia by 
means of the crossed molecular beam (CMB) technique with 
mass spectrometric (MS) detection. 

(i) O + C2H4 — > products 
(ii) O + C2H2 -¥ products 
(iii) O + CH3CCH -> products 
(iv) O + CH2CCH2 -*■ products 

In particular, the reaction mechanisms, the nature of the 
primary products of the different competing channels and 



their branching ratios have been determined (Capozza et al. 
2004; Fu et al. 2012a; Fu et al. 2012b; Leonori et al. 2012). 
From the analysis of the measured product angular and ve- 
locity distributions, the product BRs are evaluated (for the 
detailed procedure see Balucani et al. 2006 and Casavecchia 
et al. 2009). To derive the rate for each reaction channel, 
we consider the global rate constants for O + C2H2, O + 
C2H4, O + CH3CCH and O + CH2CCH2 scaled by the 
experimental BRs (see Table 1 and Table 2). For reactions 
(i) and (ii), the rate coefficients and temperature depen- 
dence are those recommended by Baulch et al. (2005) in 
the range 200-2500 K and 220-2000 K for the two reactive 
systems, respectively. For reactions (iii) and (iv), we have 
considered the rate constants and temperature dependence 
recommended by Cvetanovic (1987) in the ranges 290 - 1300 
K and 290 - 500 K, respectively. The temperatures of inter- 
est in interstellar objects are lower than those for which the 
temperature dependencies have been determined. Therefore, 
we have extrapolated the trends down to 10 K, but this is 
an approximation and important deviations are possible. 

The oxygen chemistry with hydrocarbons based on the 
new laboratory data was inserted in the UCL_CHEM chemi- 
cal model. UCL_CHEM simulates the formation of a prestel- 
lar core (Phase I) and its subsequent evolution once a star 
is born (Phase II). During Phase I gas-grain interactions 
occur due to the freeze-out of atoms and molecules on 
the grain surfaces when temperatures drop down to 10 K. 
UCL_CHEM accounts for the chemistry occurring in the gas- 
phase as well as on the grain surfaces during the free-fall 
collapse of a diffuse cloud. Species can also sublimate due to 
both non-thermal (10 K) and thermal desorption (200-300 
K) (Phase II). A more detailed description of the model can 
be found in the literature (Occhiogrosso et al. 2011). Besides 
new experimental values for the reactions listed above, the 
new part of the chemical network contains a new species, 
ketyl radical (HCCO), which is produced in the reaction O 
+ C2H2. To date, HCCO has not been detected in the in- 
terstellar medium, but its presence was already predicted 
in the past by Turner & Sears (1989) who provided an up- 
per limit for the HCCO/H2CCO ratio. We also introduce 
the vinoxy radical (CH2CHO), produced in the reaction O 
+ C2H4. Furthermore, in addition to methylacetylene, its 
isomer CH2CCH2 (allene or propadiene) is explicitly con- 
sidered. We have used either laboratory data or the rate 
coefficients included in the KIDA database (Wakelam et al. 
2012) for the reactions that involve these two hydrocarbon 
species. This distinction is not present to date in the UMIST 
database (Woodall et al. 2007) \ 

We have investigated how the revised part of the chem- 
ical model affects the trends in the hydrocarbon fractional 
abundances (relative to hydrogen nuclei) at typical densities 
in the ISM (Figure 1). Then, we have run a grid of mod- 
els representing (i) a diffuse cloud, (ii) a translucent cloud, 
(iii) a dark core and (iv) a hot core. We will analyse the 
contribution to the hydrocarbon abundances from their pri- 
mary products formed during their interactions with oxygen 
atoms. We will finally compare our theoretical hydrocarbon 



1 Just prior of submission of the present paper a new version of 
the UMIST database has been announced (McElroy et al. 2013) 
however UMIST 2012 database is not yet available. 
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Table 1. Rate constants (at 200 K, typical of a hot core) for the different competing reaction 
channels for ethylene and acetylene with oxygen atoms. The last column contains the formulae 
(taken from NIST database) adopted in order to evaluate the rate constants in column 3. 



Reaction 


BRs/% 


k/cm^mol 1 s 1 


k(T)( et /cm ti mol" I s" 1 


O + C 2 H 2 -> HCCO + H 


79 


9.9xl0~ ib 




O + C 2 H 2 ->■ CH 2 + CO 


21 


2.6xl0~ lb 


' 1.89X10- 12 (T/298) 2 - 10 c - 6535 / flT 


O + C 2 H 2 -» C 2 H + OH 











O + C 2 H 4 ->■ CH 3 



HCO 



:-!! 



l.OxlO" 13 



O + C2H4 -> CH2CHO + H 



30 



9.0xl0~ 14 



O + C 2 H 4 -> H 2 CO + CH 2 



20 



O + C 2 H 4 -> CH 2 CO + H 2 



13 



6.0X10- 14 1.01xl0- 12 (T/298) 1 ' 88 c -^5/rt 



3.9xl0~ 14 



O + C 2 H 4 -> CH3CO + H 



9.0 xlO" 1 ''' 



O + C2H4 ->■ C2H3 + OH 







Table 2. Rate constants (at 200 K, typical of a hot core) for the different competing reaction 
channels. The last column contains the formulae (taken from NIST database) adopted in order to 
evaluate the rate constants in column 3. 



Reaction 



BRs/% k/cm 3 mol 



"T7 



k(T) let /cm- i mol- 



O + CH3CCH -» C 2 H 4 + CO 



68.0 



5.2x10" 



O + CH3CCH 
O + CH3CCH 



C2H3 + HCO 



10.0 



C2H2 + H 2 CO 



10.0 



O + CH3CCH -)■ CH3CCO 



H 



5.5 



T3 - 
4.3X10" 1 " 



7.7x10" 



7.7x10" 



O + CH3CCH -> CH 3 + HCCO 



0.4 



3.1x10 



"TT7- 



O + CH2CCH2 -> C2H4 + CO 



81.5 



9.8x10" 



TT 



O + CH2CCH2 ->■ C 2 H 2 + H 2 CO 



9.0 



1.1x10" 



O + CH2CCH2 ->■ C2H3 + HCO 



7.0 



8.4x10" 



2.82x10" 



-7732/RT 



O + CH2CCH2 -> CH2CCHO + H 



1.6 



1.9x10" 



O + CH2CCH2 ->■ CH2CO + CH 2 



0.3 



3.6x10" 



abundances with those from the observations in cold and 
warm interstellar medium environments, where these species 
have been detected. 



Table 3. Loss processes for HCCO 



Reaction 



k/cm^mol 



a H + HCCO -> CO + CH 2 
a H 2 + HCCO -> products 

a O + HCCO -> CO + CO + H 



1.7xl0" 1,J 
2.8xl0~ 14 
1.6x10 



10 



2.1 The reactions O + C 2 H 2 and O + C 2 H 4 

Ethylene or ethene, C2H4, is the simplest alkene (alkenes 
are hydrocarbons with a carbon-carbon double bond), while 
acetylene or ethyne, C2H2, is the smallest alkyne (alkynes 
are hydrocarbons with a carbon-carbon triple bond). These 
two molecules have a similar structure, but their reactivity 
with oxygen atoms seems to be different (see Table 1). In 
particular, because of the larger number of atoms, ethylene 
leads to the formation of a greater number of products. Ac- 
cording to experimental results (Capozza et al. 2004, Fu et 
al. 2012a,b) in both cases the channels leading to OH for- 
mation are inefficient; this was accounted for in the UMIST 
2006 database (Woodall et al. 2007). The main channels of 
the reaction O + ethylene have been found to be those lead- 
ing to CH 3 + HCO (34%) and CH 2 CHO + H (30%) (Fu et 
al. 2012a,b). Regarding the reaction between acetylene and 
atomic oxygen, the dominant channel is the one leading to 
HCCO + H with a BR of 79% while the channel leading to 
CH2 and CO accounts for the remaining 21% (Capozza et 
al. 2004; Leonori et al. 2013). We find that HCCO plays a 
pivotal role in the trends of C2H2 abundances (see Section 
3 for more details). To date the ketyl radical has not been 
included in any database. We have also introduced a chem- 
istry network for the loss processes of HCCO as shown in 
Table 3. 



a Baulch et al. (2005) and references therein 

2.2 The reactions O + CH3CCH 

(methylacetylene) and O + CH 2 CCH 2 (allene) 

Since the experimental results by Leonori et al. (2012) and 
Balucani et al. (in preparation) have demonstrated that the 
reactions O + CH3CCH and O + CH2CCH2 are charac- 
terised by different branching ratios (and also their rate co- 
efficients are different), we have separated the two isomeric 
species in our codes. The UMIST 2006 database (Woodall 
et al. 2007) does not distinguish between the two C3H4 iso- 
mers. Therefore, it is implicitly assumed that CH3CCH and 
CH2CCH2 are formed and destroyed in the same reactions 
and with the same efficiency. Such an assumption, however, 
is not warranted. For instance, it has been suggested that 
an important contribution to the formation of both isomers 
is the reaction between methylidyne radical (CH) and ethy- 
lene (Gosavi et al. 1985; Wang et al. 1998). Nevertheless, 
recent laboratory investigations (Zhang et al. 2011) have 
shown allene to the main product with a ratio of 75-80%. 
To account separately for CH3CCH and CH2CCH2 forma- 
tion we have inserted the main path rate coefficients as rec- 
ommended by KIDA database (Wakelam et al. 2012) (see 
Table 4). In Table 2 the open channels for reactions (iii) 
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and (iv) are reported with their relative yields, as experi- 
mentally determined (Leonori et al. 2012; Balucani et al., in 
preparation) We want to point out that, at this stage, we do 
not distinguish between the two C3H3O isomers, CH3CCO 
(for the case of methylacetylene) and CH2CCHO (for the 
case of allene); however we plan to investigate their chem- 
istry in a future work. In Table 2, the main dissimilarity is 
that, in the case of allene, the channel for the formation of 
CO dominates with a ratio of 81.5%, while the BR of the re- 
action producing carbon monoxide from methylacetylene is 
only 68% (Balucani et al., in preparation). The other chan- 
nels contribute to a minor extent, especially in the case of 
the O + CH2CCH2 reaction. 



3 CHEMICAL MODELLING 

We investigated how the new reaction networks included in 
the model affect the hydrocarbon abundances by running 
simple models at typical densities for the ISM in the 10 2 - 
10 6 cm" 3 range; in particular, we simulate the case of an 
unphysical cloud at 10 K and constant density where we do 
not account for any freeze-out reaction during the collapse 
phase. We plot the fractional abundances obtained as out- 
puts from these test models for C2H2, C2H4 and the two 
C3H4 isomers as a function of the density of the cloud at 
5x10 yrs (see Figure 1). In order to highlight the changes 
in the hydrocarbon abundances, we run models with the 
same physical parameters as described above but including 
the same chemistry as in the UMIST 2006 database. In other 
words, the difference observed are only due to the new reac- 
tion scheme for reactions (i)-(iv) with all the other param- 
eters kept constant. Figure 1 shows a comparison between 
the hydrocarbon fractional abundances (with respect to the 
total number of hydrogen atoms) from the original version of 
the code (solid line) and the ones from the models containing 
the revised chemistry (dashed line). In the case of C3H4 the 
solid line refers to the output obtained from the UMIST 2006 
database that does not distinguish between the two isomers, 
while the dashed line and the dotted line are the trends for 
CH3CCH and CH2CCH2, respectively, produced from the 
updated chemistry. In the absence of freeze-out, for each se- 
lected density value we observe a discrepancy between the 
hydrocarbon fractional abundances before and after the up- 
dating, with the exception of C2H2 at 10 cm" . We empha- 
size that with this kind of comparison we are only analysing 
the effect of the new schemes on the hydrocarbon abundance 
prediction. We wish to stress that reactions (i)-(iv) are char- 
acterised by a certain activation energy and, therefore, their 
rate coefficients are very small at temperatures between 10- 
200 K. Nevertheless, they make a significant contribution 
in controlling the abundances of unsaturated hydrocarbons 
probably because of the large abundance of atomic oxygen. 
In order to understand the real influence of oxygen in 
the hydrocarbon chemistry, we modelled four different as- 
tronomical regions: a diffuse cloud, a translucent cloud, a 
dark cloud and a hot core. The following section describes 
the models in more details. 



3.1 Modelling of different astrochemical 
environments 

Table 5 lists the physical parameters used for our modelling. 
For each astronomical source we run two models: we employ 
a rate file with the same chemistry as in the UMIST 2006 
database and another one with the revised chemical net- 
work proposed in the present study. We start by considering 
the case for a diffuse cloud (Model 1), a quite extended low 
density (10 2 cm -3 ) region permeated by the ultraviolet in- 
terstellar radiation field; temperatures go up to 100 K and 
the visual extinction is around 1 mag. Because of their rel- 
atively low-visual extinction ( A„ ~ 1 mag) diffuse clouds 
are dominated by photodissociation reactions and there is 
no chemistry occurring on the grain surfaces. 

As for the case of a diffuse cloud, atoms and molecules 
do not freeze-out on the dust for a translucent cloud (Model 
2) either, an intermediate phase (in terms of physical con- 
ditions) between the diffuse and the dense medium (van 
Dishoeck 2000); the fact that they are called translucent is 
due to the presence of a bright star in their proximity that al- 
lows absorption lines to be observed. Translucent clouds are 
relatively cold regions (20-50 K) where densities can reach 
values up to 10 3 cm" 3 although there is no evidence of a 
collapsing core inside these objects. Since translucent clouds 
can be penetrated by the radiation (their A„ is between 1 
and 5 mag), their chemistry is also affected by photopro- 
cesses. Diffuse and translucent clouds are two examples of 
low-density interstellar regions. 

We also explore the case for high-density sources, such 
as dark and hot cores, where a grain-surface chemistry is 
included in order to simulate the isothermal collapse phase 
during which atoms and molecules freeze onto the surfaces of 
dust and new species can be produced. This process is more 
efficient during the formation of a hot core (a very compact 
object characterised by densities ^ 10 6 cm" 3 ) (Model 4) 
than for the case of a dark core (Model 3) where densities 
increase only up to 10 4 cm" 3 ; moreover while for the case of 
a hot core we need to model a second phase that describes all 
the heating effects that could induce molecules to thermally 
desorb (indeed temperatures rise up to 200 - 300 K), for a 
dark core we only account for the non-thermal desorption 
of species (for more details see Roberts et al. 2007) efficient 
already at 10 K. 



4 RESULTS AND DISCUSSION 

We discuss our results only for the case of a hot core (Model 
4). Results obtained as outputs from models 1, 2, 3 show 
some differences in the abundances of the species involved. 
Nevertheless we are not discussing those differences here 
since the resulting abundances are in all cases below the 
observational thresholds. Concerning the hot core model, in 
order to simulate the warm-up of icy mantles and the con- 
sequent sublimation of molecules from their surfaces, we al- 
low a time-step desorption of species as described in Viti et 
al. (2004). Based on laboratory data obtained from tem- 
perature programmed desorption experiments, Viti et al. 
(2004) categorised a list of species along with their desorp- 
tion profiles in various temperature bands. C2H2 and C2H4 
were classified as intermediate between CO-like and H2O- 
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Table 4. Paths of formation for C3H5 
database (Wakelam ct al. 2012). 



CH3CCH and CH2CCH2 based on data from KIDA 



Reaction 

H + CH3CHCH2 -► C3H5 + H 2 

C 2 H 5 + CH3CHCH2 -»■ C3H5 + C 2 H 6 

C2H3 + CH3CHCH2 -y C3H5 + C2H4 

CH 2 + CH3CHCH2 -> C3H5 + CH 3 

CH 3 + C2H3 -)■ C3H5 + H 

CCH + CH3CHCH2 -» C3H5 + C 2 H 2 



C3H5+ + NH 3 -> CH3CCH + NH 4 + 
C 2 H 6 + CH 2 CCH -t CH3CCH + C2H5 
CH 4 + CH 2 CCH -> CH3CCH + CH 3 
H 2 + CH 2 CCH -y CH3CCH + H 
C 2 H 6 + CH 2 CCH -> CH3CCH + C 2 H 5 
C2H3 + C3H5 -)■ CH3CCH + C2H4 
CH 3 + C 3 H S -y CH3CCH + CH 4 
CCH + CH3CHCH2 -> CH3CCH + C2H3 
OH + C3H5 ->■ CH3CCH + H 2 
H + C3H5 -)■ CH3CCH + H 2 
a CH + C 2 H 4 -y CH3CCH + H 



4.47X10 -13 
2.01X10" 11 
1.72 x 10- 15 
2.70xl0~ 12 
1.20xl0- 10 
2.40X10- 10 



P 



2.50 


3.50 








9.00X10" 1 " 
5.83X10- 14 
1.74xl0" 14 
1.42X10 -13 
4.30xl0~ 12 
4.00xl0~ 12 
4.03X10- 13 
2.00X10" 11 
l.OOxlO" 11 
3.00X10" 10 
7.80X10- 11 




3.30 
3.40 
2.38 




0.32 











7 



1250 


2360 
2660 









9990 

11700 

9560 

-66 


-66 











C3H5 + C3H5 -y CH2CCH2 + CH3CHCH2 1.40x10-" A^T 

C 2 H 6 + CH 2 CCH -y CH2CCH2 + C 2 H 5 5.83X10- 14 3.30 9990 

CH 4 + CH 2 CCH ->■ CH2CCH2 + CH 3 1.74X10- 14 3.40 11700 

H 2 + CH 2 CCH -y CH2CCH2 + H 1.42X10- 13 2.38 9560 

C 2 H 6 + CH 2 CCH -y CH2CCH2 + C 2 H 5 1.60X10- 12 -66 

C2H3 + C3H5 -> CH2CCH2 + C 2 H 4 2.00X10- 12 

CH 3 + C3H5 -> CH2CCH2 + CH 4 4.03X10- 13 0.32 -66 

H + C3H5 -y CH2CCH2 + H 2 3.30xl0- 10 

a CH + C2H4 -> CH2CCH2 + H 3.12X1Q- 10 



a The a value is evaluated by the product between the ratio for each channel of reaction (Zhang 
et al. 2011) and the total rate constant of 3.9xl0~ 1() cm j mor 1 s _1 

Table 5. List of physical parameters for the case of a hot-core model: visual extinction (A„), 
density (n^j ), gas temperature (T) during Phase II, efficiency of the freeze-out (fr) during Phase 
I (the percentage of mantle CO (mCO) given by the freeze-out parameter by the end of Phase I 
of the chemical model). 







A„/mag 


njj/cm 6 


T/K 


fr/% 


1. 


Diffuse cloud 


0-1 


lxlO 2 


100 





2. 


Translucent Cloud 


1-5 


lxl0 a 


20 





3. 


Dark core 


5-10 


lxlO 4 


10 


~40 


4. 


Hot core 


~500 


lxlO 7 


300 


~98 



like species. Intermediate molecules were found in the hy- 
drogenated ice layer although they were also able to escape 
via monomolecular desorption because of their moderate in- 
teraction with water ice. Since the two C3H4 isomers are 
non-polar organic species (as well as C2H2 and C2H4), we 
assumed for them a similar behaviour to the one for the 
case of acetylene and ethylene; we therefore classified methy- 
lacetylene and allene as intermediate. 

We analyse the outputs for the hydrocarbon abun- 
dances as well as the trends for the abundances of other 
species involved in the updated oxygen chemistry (see Fig- 
ure 2.). 

In particular, we focus on small patterns such as CH2, 
CH3, HCO, H2CO that have been detected as products for 
all the hydrocarbon species considered. While it is difficult 
to determine the contribution of each reaction channel to 
the abundances of all hydrocarbons, UCL_CHEM provides 
us with an analysis of the percentage of destruction and for- 
mation pathways for any selected species and we can there- 



fore draw some conclusions on the trends displayed in Figure 
2. For instance, CH2 is known as one of the reactants in the 
synthesis of hydrocarbons and its abundance is expected to 
decrease when hydrocarbons form; however we observe a rise 
in the CH2 abundances where the trends for CH3CCH and 
CH2CCH2 show a peak; this behaviour could be explained 
in the light of the fact that HCCO drops to produce CH2 
and this reaction is five orders of magnitude more efficient 
than the reaction between O and acetylene or methylacety- 
lene to form ketyl radical. Another important observation 
concerns the decrease in the ethylene abundances in order 
to produce the two C3H4 isomers, meaning that the reverse 
reaction is more likely than the channel for the formation 
of CO. CH2CO is also produced by the reaction between 
O atoms and ethylene and as a consequence its abundances 
rise when the C2H4 profile drops just after 10 6 yrs. The ma- 
jor contribution to the formation for CH3CCO is made by 
methylacetylene and, as a result, as soon as the latter species 
decreases, CH3CCO increases in abundance. Concerning all 
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Figure 1. The fractional abundances (with respect to the total number of hydrogen nuclei) for selected species as a function of the 
density. In the Top panels the solid lines represent the outputs with the same chemistry as in UMIST 2006 database and the dashed 
lines are the outputs after the updating; Bottom panel shows the trends obtained from the revised chemistry suggested in the present 
study for CH3CCH (dashed line) and CH2CCH2 (dotted line), while the solid line is the trend for both C3H4 isomers resulted from the 
previous version of the model. 



the other molecules, it is difficult to make predictions since 
they can be formed by numerous reactions, but we estimate 
a high fractional abundance for C2H3 and H2CO. This is in 
agreement with what we were expecting due to the fact that 
all hydrocarbons contribute to their formation. 

In order to validate the revised hydrocarbon chemistry, 
we evaluate the final column densities for most of the species 
involved and we compare our theoretical values with those 
from the observations towards hot cores. The column densi- 



was performed by Ridgway et al. (1976), who observed CH4 
and C2H2 towards the supergiant IRC + 10216, where a few 
years later Betz (1981) found C2H4 with a column density of 



irr 



Hydrocarbons were predicted to be abun- 



dant in the ISM but many of them were not observed at mil- 
limeter wavelengths because their symmetry forbids dipole 
rotational transitions. Some of the species listed in Table 6 
have indeed not been detected towards hot cores yet, but 
since column densities for hydrocarbons are in the 10 13 -10 15 



ties (in cm ) for a generic species Y are estimated by the cm range we expect the presence of similar amounts for 



approximation: 

N(Y) ~ X(Y) xA v x N(H 2 ), 



(1) 



where X is the fractional abundance for the species Y, A v 
is the visual extinction and N(Hv) is the column density 
of hydrogen at 1 mag and we set it equal to 1.6xl0 21 cm~ 2 . 
We also calculate the column densities obtained by the mod- 
els without the new chemistry. Results are shown in Table 
6. We observe that for all molecules, except for CH3, the 
abundances obtained as outputs from the updated version 
of the code are one order of magnitude higher than those 
evaluated from its original version. 

Hydrocarbons have been found across different astro- 
chemical environments, from planet atmospheres to dark 
clouds, in the Milky way as well as towards other Galaxies 
(Pendleton 2004). One of the first hydrocarbon detections 



C2H4 and CH2CCH2. A comparison between our theoreti- 
cal results and data taken from the observations highlights a 
general agreement within an order of magnitude; in partic- 
ular, after our updates, we achieve a greater match with the 
observed column densities of C2H2 and CH2 than if we had 
used the chemistry in the UMIST 2006 database. Concern- 
ing the two C3H3O isomers, we only calculate an upper limit 
in their column densities since we do not account for a com- 
plete reaction network for their loss mechanisms. We would 
also like to stress that we have only qualitatively modelled 
a hot core without looking at a specific high-mass source 
therefore our theoretical column densities could be refined. 
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Figure 2. The fractional abundances (with respect to the total number of hydrogen nuclei) for selected species as a function of time. 



Table 6. Comparison between observational and theoretical column densities (in cm -2 ) for se- 
lected species. Starting from the left, the first column shows results from the original version of 
the code; on the second column we report the molecular column densities calculated after our 
updating. 



Molecule 


UMIST 2006 


Present Work 


Observations 


Source 


References 


C 2 H 2 


2.70xl0 lb 


4.08X10 14 


5.50X10 14 


SgrA 


Lacy ct al. (1989) 


C 2 H 4 


1.93X10 13 


1.26X10 14 


- 


- 


- 


CH3CCH 


l.llxlO 15 


6.18X10 15 


1.20X10 15 


Orion-KL 


Wang et al. (1993) 


CH2CCH2 


not present 


1.78X10 16 


- 


- 


- 


CH 2 


2.09X10 15 


4.88X10 14 


6.60X10 13 


Orion-KL 


Hollis et al. (1995) 


CH 3 


2.70X10 15 


7.73X10 13 


8.00X10 14 


SgrA 


Feuchtgruber et al. (2000) 


HCCO 


not present 


4.68x10 s 


- 


- 


- 


CH 2 CO 


4.25X10 14 


2.55X10 15 


l.OOxlO 14 


SgrB2 


Turner et al. (1977) 


HCO 


5.62X10 12 


5.13X10 12 


1.03X10 12 


W3 


Snyder et al. (1976) 


H 2 CO 


6.55X10 14 


1.85X10 14 


8.32X10 14 


Orion-KL 


Wright et al. (1996) 


CH3CO 


not present 


9.13X10 14 


- 


- 


- 


CH3CCO 


not present 


5.80 xlO 15 


- 


- 


- 


C2H3 


5.44X10 13 


1.30X10 13 


- 


- 


- 



5 CONCLUSIONS 

In this paper we have presented how a more accurate chemi- 
cal network representing the reactions of atomic oxygen with 
several unsaturated hydrocarbons affects the abundances of 
simple organic species in ISM. We mostly focus on the mod- 
elling of hot ISM regions where these variations were found 
to be more evident; however colder regions such as diffuse 
and translucent clouds or a dark core have also been inves- 
tigated. A first conclusion is that towards hot cores/corinos 



atomic oxygen easily degrades unsaturated hydrocarbons di- 
rectly to CO or to its precursor species (such as HCCO or 
HCO) and destroys the double or triple bond of alkenes and 
alkynes. Therefore, environments rich in atomic oxygen at a 
relatively high temperature are not expected to be charac- 
terised by the presence of large unsaturated hydrocarbons 
or polycyclic aromatic hydrocarbons. On the contrary, in O- 
poor and C-rich objects, hydrocarbon growth can occur to 
a large extent. In colder objects, the effect of the reactions 
of atomic oxygen with small unsaturated hydrocarbons is 
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reduced because of the presence of a small activation en- 
ergy for these reactions. Moreover we highlight the main 
contribution from other reaction channels that are not yet 
included in the databases available: ketyl and vinoxy rad- 
icals are among the new significant species that should be 
considered in the modelling. New attempts to detect them 
should be pursued and a chemical scheme for the vinoxy rad- 
ical loss pathways devised. Furthermore, we have established 
that structural isomers, such as methylacetylene and allene, 
should be treated separately when they manifest a different 
chemical behaviour. A more general conclusion is that, if wc 
wish to fully understand the formation of the numerous or- 
ganic species detected in ISM, the intricate chemistry that 
leads to their formation should be treated in detail. 

We finally plan to refine our theoretical results by mod- 
elling specific sources where hydrocarbons have been found. 
In the present study we calculated the hydrocarbons column 
densities that match those from the observations within an 
order of magnitude although we were unable to make a com- 
parison for all the species considered: most hydrocarbons are 
indeed non-polar molecules and therefore they cannot be de- 
tected. Despite this constraint, we aim to estimate the col- 
umn densities for the case of symmetric hydrocarbons based 
on the abundances of their derivative species; for instance, 
as already stated by Quan & Herbst (2007), the amount of 
allene in a source could be related to the abundances of its 
derivative cyanoallene that, contrary to allene, is detectable. 
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